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ABSTRACT 

We suggest that the observed long-period Alfven waves in the solar wind may be 
generated in the solar interior due to the pulsation of the Sun in the fundamental 
radial mode. The period of this pulsation is about 1 hour. The pulsation causes a 
periodical variation of density and large-scale magnetic field, this affecting the Alfven 
speed in the solar interior. Consequently the Alfven waves with the half frequency of 
pulsation (i.e. with the double period) can be parametrically amplified in the interior 
below the convection zone due to the recently suggested swing wave- wave interaction. 
Therefore the amplified Alfven waves have periods of several hours. The waves can 
propagate upwards through the convection zone to the solar atmosphere and cause 
the observed long-period Alfven oscillations in the solar wind. 
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1 INTRODUCTION 

Recent in situ observations of solar polar regions by Ulysses 
spacecraft have shown the presence of long-period (peri- 
ods of several hou r s) outwardly propa gating Alfven waves 
', JSmith et al.lll995t iBaloeh et ai]ll995l) . The Alfven waves 
• may play a significant role in solar wind acceleration, there- 
fore it is very important to understand the process of their 
generation. As the observed Alfven waves show predomi- 
nantly outward propagation, they probably are of solar ori- 
gin. Unfortunately no clear physical mechanism explaining 
the generation of such long-period Alfven waves has been 
suggested so far. 

It is an interesting fact, that the low- frequency os- 
cillations have been observed also in solar lower at- 
mosphere iMer kulenko et al.lll 988l:lMerkulenko fc Polvakovl 
ll99ll:lMashnich fc Bashkirtsevlll999l) . The Amplitude of ob- 
served oscillations grows from the center of the solar disc 
to the limb i.e. the waves are rather tangential than radial. 
Therefore the authors suggested that they can be Alfven 
waves propagating upwards from the solar interior. Thus 
the energy source of these Alfven waves can be located well 
below the solar surface. 

Here we suppose that the long-period Alfven waves can 
be generated in the solar interior by pulsation either in 
the fundamental radial mode and/or in low frequency g- 
modes through the recently suggested swing wave-wave in- 



teraction dZaaarashvilil 


20011: IZaaarashvili & Robert J2002a: 


IZaaarashvili et alJl2002 


) . The period of the solar fundamen- 



tal radial mode is ~ 1 hour; periods of g-modes can be 



* E-mail: temury@genao.org 
f E-mail: gbelvedere@ct.astro.it 



longer. The presence of a large-scale weak magnetic field in 
the interior can not affect the pulsation significantly. How- 
ever the pulsation causes the periodical variation of some 
medium parameters (such as density and magnetic field) and 
consequently of the Alfven speed in the interior. But the pe- 
riodical variation o f Alfven speed either due to soun d waves 
JZaaarashvilil 120011 IZaaarashvili fc Roberts! l2002bl) or fast 
magnetosonic waves JZaaarashvili fc Robertsn 2002a) leads 
to the amplification of Alfven waves with twice the period of 
compressible waves. By this mechanism, it has been shown 
that the pulsation in the fundamental mode of solar-like 
stars in binary systems, which is excited by the compan- 
ion's gravity, may lead to the amplificat ion of double period 
Alfven waves iZaaarashvili et alJ l2002). Thus the question 
naturally arises whether a similar process may occur in the 
Sun. In this case, the solar fundamental pulsation may am- 
plify torsional Alfven waves in the interior with the double 
period i.e. ~ 2 hours, while g-modes may amplify Alfven 
waves with longer periods. 

The generated Alfven waves may propagate up- 
wards and cause the observed oscillations in the photo- 
sphere/chromosphere and the solar wind. It is a very in- 
triguing fact that one of ob served oscillation period in the 
chromosphere is ~ 2 hours jMerkulenko et allll988l) . 

In this letter we show how the solar fundamental radial 
pulsation may amplify the torsional Alfven waves in the in- 
terior. As to the coupling of g-modes and Alfven waves, it 
is more complicate to describe it mathematically, therefore 
we will study it in future. 
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2 STATEMENT OF THE PROBLEM AND 
DEVELOPMENTS 

We start from the ideal magnetohydrodynamic (MHD) 
equations: 

^ + V(pv)=0, (1) 

dv 1 
p— + p(v-V)v = -Vp + — (VxB)xB- P V0, (2) 

f =Vx(vxB), (3) 
V 2 = A-wGp, (4) 

P = Po(g)\ (5) 

where p is the medium density, p is the pressure, v is the 
velocity, B is the magnetic field, G is the gravitational con- 
stant, (f> is the gravitational potential and 7 is the ratio of 
specific heats. We consider a medium with zero viscosity, in- 
finite conductivity and negligible displacement current. We 
argue that the coupling between pulsation and Alfven waves 
occurs below the convection zone, therefore we do not con- 
sider convective effects and consequently the Q-effect term 
in the Mean Field Electrodynamics induction equation (see 
e.g. Belvedere f 985) which represents a mean electromotive 
force generated by cyclonic convection. Here and in the re- 
maining part of the paper we also neglect the rotational 
effects. 

2.1 Fundamental pulsation of the Sun 

The Sun is considered to be in hydrostatic equilibrium so 
that its own gravity is balanced by the pressure gradient. 

A slight deviation from the equilibrium leads to an oscil- 
lation which may be studied by the linear perturbation the- 
ory. Stellar oscillations (in the absence of magnetic field) can 
be divided into two classes: high-frequency p-modes or pres- 
sure modes, whose restoring force is the pressure gradient, 
and low-frequency g-modes, where buoyancy is the restor- 
ing force. There are also interm ediate oscillations sometimes 
called f-modes JCowlindlfQllh . 

The boundedness of any system usually leads to oscil- 
lation in the first eigenfrequency which is also called the 
fundamental frequency of the system. The mechanical ana- 
log of the fundamental oscillation is the tuning fork in the 
case of an impulsive force and the musical trumpet in the 
case of a continuous force. The oscillation of a star in the 
fundamental mode yields a wavelength comparable with the 
stellar radius. Therefore the fundamental frequency evalu- 
ated for an homoge neous self-gravitating sphere will be of 
order of y/GM/Rt^Cs/R, where R and M are the radius 
and the mass of star and c s — -y/TPo/po is the mean sound 
speed in the interior. However the real equilibrium density is 
not uniform throughout the Sun, therefore the fundamental 
period is shorter ~ f hour. Unfortunately no clear obser- 
vational evidence of the fundamental pulsation was found, 
which is probably due to the acoustic cutoff in the stratified 
atmosphere. The acoustic cutoff frequency in an isothermal 
atmosphere is c s /2Ao, where Ao is pressure scale height. Un- 
der conditions of the photosphere/chromosphere Ao « 125 
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Figure 1. Schematic picture of generation and propagation of low 
frequency Alfven waves in the Sun. The Alfven waves are gener- 
ated in the solar interior due to the pulsation of the Sun in the 
fundamental radial mode. Then they propagate along magnetic 
field lines towards the corona. 

km. Then for the typical photospheric sound speed c s = 7.5 
km/s the cutoff frequency is 0.03 s , which gives the cutoff 
period of 210 s. This means that the sound waves with longer 
periods, such as the fundamental radial pressure mode, can 
not pass through the photosphere. Probably this is a reason 
why the low-frequency pressure modes are not seen in oscil- 
lation spectr um (but see a possib le evidence of low frequency 
pulsation in ferown et al.lll97gli '). 

As spherical coordinates lead to a more complicated 
mathematical formalism, we adopt a cylindrical coordinate 
system (r, cf>, z) and for simplicity consider only the axisym- 
metric problem so everywhere d/d<j> = is assumed. We 
suppose that an uniform poloidal magnetic field, Bo, exists 
inside the star directed along the z axis (see Fig.l). The 
large-scale magnetic field is generated either by classical dy- 
namo or other mechanisms, but this is beyond the scope of 
the present paper. We just assume that the magnetic field 
exists. This field is considered to be relatively weak, so it 
does not affect the fundamental pulsation significantly. 

In the adopted cylindrical coordinate system, we ap- 
proximate the fundamental solar oscillation as a radial pul- 
sation of the cylinder, thus neglect the motion along the z 
axis. Then the linearised equations (l)-(5) can be split into 
radial and azimuthal components, where the radial compo- 
nent corresponds to pulsations and the azimuthal compo- 
nent to torsional Alfven waves. Torsional Alfven waves have 
azimuthal velocity polarisation and propagate along the 2 
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axis. As we will see later, the wavelength of the Alfven waves 
is much smaller as compared to the solar radius (due to small 
Alfven speed in the interior) , therefere the z velocity compo- 
nent of the global pulsation, which has a spatial scale compa- 
rable to the solar radius, can be considered as homogeneous 
in Alfven spatial scale and therefore cannot significantly af- 
fect the Alfven wave dynamics. Therefore, the motion along 
the z axis due to the global pulsation may only complicate 
the calculations and thus is neglected as mentioned above. 

The radial part of linearised equations (l)-(5) in cylin- 
drical coordinates is 
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while the azimuthal part is 

db<f, _ du<t, 

dt dz 
du^ _ Bp db$ 

dt Aitpo dz ' 



(10) 
(11) 



where u r , u<f,, b z , b$, 8p and 5p are the velocity, magnetic 
field, density and pressure perturbations respectively. Equa- 
tions (6)- (9) govern the linear radial pulsation of the star, 
while equations (10)-(11) describe the linear torsional Alfven 
waves. 

In general, stellar radial adiabatic pulsations can be rep- 
resented as standing spherical waves l|Coxlll9 8fO with a lin- 
ear radial velocity field 



u r = aF(r)sin(ui n t), 



(12) 



where uj n is the eigenfrequency, F{r) is the eigenfunction 
and a is the pulsation amplitude. The expression of the 
eigenfunction F(r) depends on the spatial profiles of the 
unperturbed physical quantities and may be represented by 
some combination of spherical Bessel functions. However, 
here we are not interested in specific pulsation functional 
forms, therefore we retain its general form F(r), which then 
may be specified by choosing the spatial distribution of the 
unperturbed quantities throughout the Sun. 

Then, using expression (12), the linearised continuity 
and induction equations give the density and magnetic field 
perturbations as: 



Sp = — cos(uj n t)Fp(r), 

Oi 

b z = — cos(u> n t)F b (r), 
where 

F P (r) = Po(r)^r^ + 
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dr 
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(13) 
(14) 

(15) 
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where po(r) is the unperturbed density. 

Equations (13)-(14) show that the radial pulsation leads 



to local periodical variation of density and magnetic field 
throughout the Sun, which causes the periodical variation 
of the Alfven speed at each level r. In next section we show 
that it determines an exponential amplification of torsional 
Alfven waves with half the frequency of pulsation. 

2.2 Resonant torsional Alfven waves in the solar 
interior 

The radial pulsation affects the torsional Alfven waves 
through nonlinear interaction and then equations (10)-(11) 
can be rewritten as 



db$ _ ^ du^ _ d(u r bj,) 

dt z dz dr ' 

duj, ^ u r d{ru4>) _ (Bo + b z ) dbj, 
dt r dr 47r(po + Sp) dz 



(17) 
(18) 



Thus equations (17)- (18) describe the evolution of 
Alfven waves in presence of the global pulsation. Since the 
energy associated to this pulsation is clearly much larger 
than the energy stored in Alfven waves, the back reaction 
of Alfven waves can be ignored (at least at the initial phase 
of the wave evolution), which means that the amplitude of 
the pulsation remains constant. 

Equations (17)-(18) can be further simplified by taking 
the Alfven oscillations in narrow cylindrical shells, then the 
unperturbed density po and magnetic field -Bo can be con- 
sidered as homogeneous (thus r will stand as a parameter in 
the equations) and equations (13)- (14) give 

Sp _ b z 
Po B ' 

Then equations (17)-(18) lead to the second order dif- 
ferential equation 
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Using the Fourier transform of u$ with z-dependence 

(21) 

equation (20) is rewritten as 
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= 0. (22) 
Now, using the well known substitution of function 

_1 r 8(r. r ) , 

ut = u 4 ,e 2 J ~^- d \ (23) 

the term with the first time derivative can be removed and 
with expressions (12)-(14) we finally get 
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Without the pulsation (i.e. taking a=0), equation (24) 
describes the usual dispersion relation of Alfven waves. 
But, in presence of the forcing introduced by the pulsa- 
tion (cv^O), equation (24) is the well known Mathieu equa- 
tion. The solution of this equation wi th frequency o) n /2 
has a n exponentially growing character flLandau fc L ifshitz 
ll988|), thus the main resonant solution occurs w hen 
JZaaarashvili fc Roberti]l2002at IZaaarashvili et aljEoolj) 



(26) 



where uia is the frequency of torsional Alfven waves. Under 
condition (26) the solution of equation (24) is 



u$(t) = uoe~ 



cos — t=Psm — 
2 ^ 2 



(27) 
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where uo = u(0) and the phase sign depends on a^I'(r); it is 
+ for negative a^(r) and — for positive a^(r). Note that 
the solution has a resonant character within the frequency 
interval 



(28) 



Thus the radial pulsation leads to an exponential am- 
plification of torsional Alfven waves with the half frequency 
hu> n - The growth rate depends on the amplitude a and spa- 
tial structure of the pulsation eigenfunction F(r), which in 
turn depends on the radial structure of density po and mag- 
netic field Bo. The resonant condition (26) imposes a restric- 
tion on the wave number k z . In other words, the pulsation 
"picks up" the harmonics of torsional Alfven waves with a 
certain wavelength at each radial distance r. 

Numerical simulation of equations (17)-(18) also shows 
the resonant amplification of torsional Alfven waves. For 
simplicity we take the function F in expression (12) as con- 
stant and look for the temporal evolution of the spatial 
Fourier harmonics of &</, and u^. The results are shown on 
Fig. 2. The exponential amplification of the Alfven wave spa- 
tial Fourier harmonics (upper plots) with the double period 
of pulsation is clearly seen. 

As the main resonance occurs at the half frequency of 
pulsation (see equation (26)), the resonant Alfven waves in 
the solar interior will have a period 



T A = 2T ~2 hours, 



(29) 



where To ~ 1 hour is the pulsation period. 

Due to the resonant range (28), the generated Alfven 
waves will have not a fixed period, but many periods in 
that range. The width of range depends on the amplitude 
of pulsation. Thus the period of generated Alfven waves in 
the solar interior may be in the range of 



Ta ~ 1.5 — 3 hours. 



(30) 



From expression (26) we may calculate the wavelength 
of the generated Alfven waves. Using the dispersion rela- 
tion of the sound waves, we have \/R ~ va/c s , where A 
is the wavelength of the resonant Alfven waves and va = 
Bo I V4ttPo is the mean Alfven speed in the interior. As the 
plasma /3 (/3 « cI/v'a) ls higher in the solar interior, the 
wavelength of the Alfven waves will be shorter comparing 
to the solar radius. Taking c s /va = 100, we get a value of 
this wavelength in the solar interior such as A ~ 7 000 km. 

It must be mentioned that the amplification of the 
Alfven waves due to the fundamental pulsation can occur in 
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Figure 2. Resonant amplification of Alfven waves due to the 
radial pulsation. The plot below is the velocity of the fundamental 
pulsation. The plots above are the magnetic field and the velocity 
components of the torsional Alfven waves. These waves have the 
double period of pulsation. The energy of pulsation is considered 
to be much larger than the energy of the Alfven waves, therefore 
the amplitude of pulsation remains constant. 



the solar interior below the convection zone, as the motions 
in the convection zone may destroy the resonance process. 

The amplified Alfven waves can propagate along the 
magnetic field lines through the convection zone and reach 
the solar atmosphere (see the schematic picture in Fig. 1). 
The configuration of large-scale magnetic field is very sim- 
plified as the complication due to the convection is com- 
pletely ignored. In fact the stochastic small scale magnetic 
fields in the convection zone may make the wave propaga- 
tion difficult. But the magnetic fields in sunspots and in 
regions between large convective cells are predominantly 
large scale and they may easily guide the waves through 
the convection zone. Therefore the Alfven waves probably 
can be observed in sunspot latitudes and in polar regions 
where the dipole component of magnetic field dominates. 
It is an interesting fact that observed low-frequency os- 
cillations found in photosphere, chromos phere and promi- 
nence by iMashnich fc Bashkirtsevl l)l999f) have the period 
of ~ 1.75 hours at sunspot (25°-30°) and higher (60°) lati- 
tudes. The observed oscillations are rather transversal than 
the radial and can be associated to the Alfven waves prop- 
agating from the solar interior. Also long-period transver- 
sal oscillations with the period of ~ 2 hours have been 
observed in chromospheric spectral lines jMerkulenko et all 
ll988UMerkulenko fc Polvakovll99Jh and again can be inter- 
preted as Alfven waves. The lo ng-period oscillations recently 
observed in solar prominence jTerradas et aUl2002ft can be 
also due to the Alfven waves propagating from below. It is 
also intriguing fact that long-period Alfven waves are ob- 
served in fast solar wind at high latit udes where dipolar 
magnetic field com ponent is dominant iSmith et alJll99f3 : 
Balo gh et al"lll995F) . It is quiet possible that photospheric 
and chromospheric Alfven waves may propagate into the 
corona and solar wind. The amplitude of long-period oscil- 
lations jMashnich fc Bashkirtsevl [l999'l in the photosphere 
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is ~ 0.1 km/s and in the chromosphere ~ 0.25 km/s. If these 
oscillations propagate into the low density corona then the 
amplitudes must be increased as the wave energy density re- 
mains approximately constant. The wave energy density is of 
order ~ pit 2 , where p is the background density and u is the 
wave velocity. The particle density in the photosphere is 10 s 
times more than in the corona. Thus in order to retain the 
constant energy density, the wave velocity must increase 10 4 
times i.e. to the value ~ 10 3 km/s. Which means that the 
waves become intrinsically non-linear as the Alfven speed 
has the similar value in the corona. It is the fact that solar 
wind Alfven waves are nonlinear 5V /V ~ 8B/B ~ 1. 

Another question is the energy balance between solar 
global modes and observed Alfven waves. The energy stored 
in the fundamental pulsation is huge as the whole Sun takes 
part in the oscillation. Also the ratio between the densities 
in the interior (say the base of convection zone) and the pho- 
tosphere is ~ 10 7 . Therefore even low amplitude pulsation 
can be responsible for observed Alfven wave energy. 

On another hand, as we mentioned above there is 
no clear observatio nal evidence of solar global pulsation. 
iBrown et alJ l)l978h reported the variation of solar diameter 
with the periods around 1 hour. But since the pulsation was 
not clearly observed in solar oscillation spectrum. There can 
be two main reasons why the fundamental pressure mode is 
absent at the photosphere. Either it has very low amplitude 
or it can not pass the photosphere due to the stratification 
cutoff (see above). On another hand, the continuous energy 
transfer from the pulsation into the Alfven waves in the in- 
terior may significantly lower the amplitude of pulsation at 
the surface. 

As the mean Alfven speed varies from the solar interior 
to the corona, the wavelength and the period of these up- 
wards propagating Alfven waves will change along their way 
to the corona. Therefore the observed Alfven oscillations in 
the solar wind will have a longer wavelength than the Alfven 
waves in the interior. 

It must be mentioned that many simplifications (such 
as the cylindrical symmetry instead of the spherical one, the 
spatially homogeneous distribution of physical quantities) 
have been made to outline the physics of the Alfven wave 
generation in a simple way. Of course, the real situation is 
much more complicated, but the physical meaning of the 
mechanism will remain the same. Future exstensive numer- 
ical simulation is needed in order to take into account the 
real distribution of physical quantities, but this is beyond 
the scope of this letter. 



3 CONCLUSIONS 

We suggest that outwardly propagating long-period Alfven 
waves, which have been observed in solar po lar regions by 
Ulysses iSmith et alll993lBaJogh et alll995ft can be gener- 
ated in the solar interior due to the solar fundamental radial 
pulsation. It is shown that this radial pulsation parametri- 
cally amplifies torsional Alfven waves. The main resonance 
occurs for Alfven waves with the double period of pulsation, 
although the process has a resonant range of frequencies, 
whose width depends on the amplitude of pulsation. There- 
fore the solar fundamental radial pulsation, with a period 
of ~ 1 hours, will amplify Alfven waves with a period of ~ 



1.5-3 hours. These amplified Alfven waves may propagate 
upwards to the solar corona and cause the observed mag- 
netic field and velocity oscillations in the solar wind. 

It must be mentioned that low-frequency g-modes may 
also amplify Alfven waves with longer periods, but we do not 
consider this process in this letter. However it is really worth 
to study the process in future, as the existence of g-modes 
in the solar interior could be revealed through Alfven oscil- 
lations in the photosphere/chromosphere and solar wind. 
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